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Abstract 

Rare kaon processes appear to be particularly suitable to study the extensions 
of the standard model, especially if the possibility for eventual direct evidence 
becomes unlikely. In this review, we discuss processes that are important as a 
test of either the standard model or supergravity. Moreover, some of these are 
important even for both the standard model and for supergravity. 

Particular attention is paid to the reduction of uncertainties in the calculation, 
especially the ones coming from the confinement effects. Recent approaches, such 
as chiral perturbation theory, the large A'^c-expansion, QCD sum rules and lattice 
QCD, are discussed. This is found to be the best strategy in view of the fact 
that supersymmetric effects are rather tiny. 
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1 Standard Model and Super symmetry 



1.1 Introduction 

Three basic interactions (electromagnetic, weak, and strong) seem to be reasonably (even 
very well) described by a theory known as the standard moc?e/|jll-|jTO|- The fourth known 
interaction, gravitation, is not included. Its classical version, the general theory of relativity, 
describes macroscopic phenomena rather well; the quantized theory of gravity, in spite of 
many attempts, still does not exist as a self-consistent theory. Fortunately, gravitational 
interactions at energies of present experiments do not play any important role. Therefore, 
the standard model may be tested to a high degree of accuracy. 

Up to now there has been no confirmed disagreement between the standard model and 
experiment. The standard model describes many phenomena very well; actually, the de- 
scription is too good for the taste of many theoretical physicists. This 'unhappiness' arises 
from the answer to our question about the standard model: the standard model is not a 
satisfactory theory for the reasons to be discussed later. 

Theoretically, one tries to construct theories which go beyond the standard model and 
improve the drawbacks of the standard model in some aspects. This includes a desirable 
reducing of the number of parameters, improvement in convergence properties and, particu- 
larly, the understanding of the Higgs sector which, in itself, at least in the minimal version 
of the standard model, shows unpleasant behavior at the quantum level. Clearly, theories 
which go beyond the standard model [0, bring new physics into play - usually a plethora 
of new particles. Naively, this looks very promising because new particles may, at least in 
principle, be detected in one or another way. However, almost perfect agreement between 
the standard model and experiment requires that new physics should have very little impact 
on low-energy phenomenology. This means that new particles are rather heavy. Therefore, 
their direct production and detection at large accelerators might finally appear very difficult, 
especially if the present upper limits on masses of new particles become considerably larger. 

An alternative but complementary way of searching for new physics are low-energy ex- 
periments with a high degree of accuracy. This way is based on the following points. 

(i) There are processes that appear at higher order in electroweak interactions and which 
are very sensitive to 'impurities' caused by new physics. A special role is here played by 
rare kaon processes. These include rare kaon decays, CP violation in the kaon system, e'/e 
parameter, etc. 

(ii) Significant progress has been made in reducing the uncertainties in theoretical pre- 
dictions in the standard model. These uncertainties are usually connected with the treating 
of hadrons, i.e., with the problem of QCD confinement. A few new, rather sophisticated ap- 
proaches have been addressed to this problem: QCD lattice calculations |jl2|, the large- iVc 
approach ||14||, and the QCD-duality approach ||15||. 



Since the effects of new physics appear to be rather tiny, the importance of reducing the 
uncertainties in calculations and/or errors in experiments is conditio sine qua non for the 
future progress. 

A number of experiments are presently running at BNL, KEK, CERN and FNAL. Kaon 
beams are usually produced from high-energy proton beams colliding on fixed targets. Ex- 
ception is the CP-LEAR project at CERN where low-energy pp collisions are used to study 
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kaons in the final state, especially CP violation and CPT tests. 

A clean, intense source of K-mesons is expected in 0-factories, like DA$NE in Frascati, 
which is to be built in the next few years [0]. 

A few proposals for K-factories have been considered in the last few years, like Euro- 
pean Hadron Facility Jl^ (abbandoned) and TRIUMPH [|l^] with extremely intense beams 
(planned after 1995). 



Rare Kaon Decays 

The kaon system has proved to be the graveyard of many wrong theories. The history 
goes from the t — 9 puzzle, CP violation, the A J =1/2 rule, the absence of AS" = neutral 
currents to the present search for lepton-flavor violation and the existence of the fourth 
generation. Last but not least, new ideas, e.g., composite models and/or supersymmetry 
have to pass hard tests in kaon physics. 

The standard model (SM), based on the gauged SU{?>) x SU(2) x U{1) theory with the 
minimal Higgs sector, can be proved/disproved in the following ways: 

(i) Going to higher energies and looking whether the SM still works; this is basically the 
approach of collider physics. 

(ii) Precision measurements: measuring the parameters of the SM highly accurate and 
looking for discrepancies. 

(iii) Searching for processes which are suppressed or forbidden in the standard model. 
This basically defines the physics of rare processes, which are the subject of this paper. 

An example of forbidden process is 

KL^lie, (1.1) 
and, clearly, its observation at any level would be a clean signal for new physics. The present 



experimental limit is [[19|j 

5i?(is:L /i e) < 0.94 X 10"^°. (1.2) 

Other examples of K decays that are forbidden in the SM, but are predicted in new models 
are 

7r+/i+e-, 
Kl 7r°/ie, 

K+ TT+XO, (1.3) 

where is a light scalar or pseudoscalar. 

There are many processes in kaon physics that are strongly suppressed in the standard 
model, but are not totally forbidden. We list typical examples. 

The first process, Kl — > e"'"e~, is especially sensitive to the Higgs sector of new physics 
since, unlike in the SM, there would be no large helicity suppression. 

The Kl — >■ 7r°e"'"e~ process is suppressed by a factor ~ 10~^ owing to the GIM mechanism 
and, in addition, by a factor of ~ 10~^ owing to CP conservation. 
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The Tx^uu is probably the cleanest test of higher-order electroweak interactions in 

the standard model. This particular process belongs to the class of more general K decays 
of the type 

ir+ ^ 7r+ + nothing, (1.4) 

where 'nothing' denotes any light neutral particle(s) that cannot be detected. 
Particularly interesting is the process 

K+ 7r++?+?'. (1.5) 

This channel receives contributions from standard decays (i.e., — > tx^vv ) and from the 
new physics beyond the standard model. Being practically a short-distance phenomenon, 
it has the advantage with respect to processes such as Ki — > or 7r+e+e~, where 

long-distance contributions enter the game. Also, any mechanism for lepton flavor violation 
(horizontal gauge bosons, leptoquarks, etc.) would lead to — > vr+z/eZ/^. In addition, flavor- 
changing neutral currents in N=l supergravity would give a contribution to — > tt^vv in 
addition to the SM contribution. 

Finally, the direct evidence for superparticles would be possible in processes such as 

K+ ^ 7r+77, (1.6) 

provided the photinos 7 are light enough to be produced. Last but not least, the existence 
of the fourth generation would certainly influence this rare decay. 



1.2 Limits of the Standard Model and Supersymmetry 

The standard model is a theory with a lot of parameters which are not explained by the 
model itself. These include lepton and quark masses, Yukawa couplings, mixing angles, etc., 
whose pattern is not understood. Why have left-handed fermions so far appeared in SU (2) 
doublets and right-handed ones in SU{2) singlets? Why are there three colors, and why is 
the electric charge quantized? The family problem is not understood: why three generations 
of quarks and leptons?, etc. 

Searching for answers one is forced to go beyond the standard model. Among different 
theories, supersymmetry is particularly interesting||20||-[^ ; it is the symmetry of bosons and 
fermions, whose theoretical motivation we discuss in the following. 

(i) Experience teaches us that nature obviously prefers gauge theories and supersymmetry 
is the next logical gauge theory. 

(ii) The spin degree of freedom is naturally contained only in supersymmetry. 

(iii) Supersymmetric theories are much better behaved mathematically, and are even 
flnite in some cases. 

(iv) Locally gauged supersymmetry relates the SUSY generators to the generators of 
space-time transformations leading to supergravity, i.e., to a natural coupling of SUSY to 
gravity. 

The standard model can hardly be considered as a fundamental theory since it contains, 
e. g., QED, which is not an asymptotically free theory. Its interactions must become strong 
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at some higher energy scale. This is considered as a hint that one has to treat the standard 
model rather as a low-energy effective theory of a more fundamental one. 



The Higgs sector 

The Higgs sector gives additional hints which cast doubt upon the fundamental character 
of the standard model. The minimal Higgs sector needed to construct the standard model 
contains one complex Higgs doublet $ with the potential 

= A[$t$ - (1.7) 
After spontaneous symmetry breaking, this leads to one physical Higgs field, with the mass 

= ^A, (1.8) 

with an unknown coupling constant A. 

As far as the Higgs mass is concerned, it cannot be taken arbitrarily large. Radiative 
corrections lead to a lower bound: 

Mh > 6.6 GeV. (1.9) 
The various analyses |l], ^ show that the upper limit on the Higgs mass is 

Mh < 1 TeV. (1.10) 

Once the upper limit, Mh < a few TeV, has been set, the huge radiative corrections to the 
mass of the Higgs boson provide an additional hint that one needs a new physics beyond 
the standard model. These corrections arise, e.g., from different loops, which turn out to 
be quadratically divergent, as is typical of scalar theory, i.e., 6Mfj ~ A^. Since it would be 
unnatural to expect the corrections to be larger than the upper limit on the Higgs mass, one 
expects the new physics to give an effective cutoff scale below a few TeV . 

The idea is to use a higher symmetry to eliminate quadratic divergences in the Higgs 
mass corrections. Here supersymmetry appears to be the proper theory because it can in 
principle eliminate the problem of quadratic divergences. The mechanism is rather simple: 
In supersymmetry the loops of normal particles are always accompanied by the loops of 
superpartners. The extra minus sign, appearing because of the fermion loop, leads basically 
to the cancellation of divergences. This nice property persists as long as the imposed super- 
symmetry is exact. Since supersymmetry necessarily has to be broken, the requirement of 
approximate cancellation of divergences imposes a constraint on the masses of superparticles. 

The Higgs mechanism brings in the vacuum expectation value v which determines Mw 
and fermion masses: it is u ~ 250 GeV . This value is not derived in the standard model 
itself and needs a fundamental explanation. Any attempt made so far to explain v in terms 
of higher symmetries brings us to the fine-tuning and the hierarchy problem. 

Let us call Hweak a scale at which SU{2) x U{1) breaking takes place, and M a scale 
where a fundamental theory becomes relevant. Typically, this is a very high scale, e.g.. 
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Mgut ~ 10"^"^ GeV, Mpianck ~ lO"*^^ GeV , etc. Now, in the fundamental theory, the mass of 
the Higgs particle is evaluated at the fundamental scale M, and M^(M) has to be scaled 
down to the weak scale [iweak- Typically, the expression is 



I'M 

Ml{yi^,ak) = Ml{M) + const, / dv^ + ^{M)g^ + (!)(/), (1.11) 



where 3fJ(M) grows almost as InM when M —>■ oo, and g is a coupling constant. The second 
term in ( |1 . 1 1| ) diverges quadratically as M — > oo. 

In order to have M^(/x^eaA;) M"^, one has to fine-tune the parameter M'jj{M) extremely 
accurately to cancel the second term in ( |1 . 1 1| ) which is if order M^. This is known as the 
fine-tuning problem or the naturalness problem. Suppose, for example, that the standard 
model is valid up to a GUT scale of 10^^ GeV or even up to the Planck scale of 10^^ GeV . If 
these two scales plus the weak scale were input into the theory, the Higgs mass would have 
to be chosen with an accuracy of 10~^^ compared with Mpianck- Clearly, the natural value 
for M]j{fi^eak) is ~ 0{M^). Therefore, the fundamental question is actually the hierarchy 
problem: why is fiweak |^, |T^? Even if the question of fine-tuning were solved in a 

satisfactory way, the hierarchy problem would still have to be solved. 

As we have already mentioned, one has to impose constraints on possible supersymmetry 
breaking to preserve the cancellation of divergences. Such a symmetry breaking is called soft. 
An example of soft symmetry breakdown is a spontaneous symmetry breaking. Clearly, with 
softly broken SUSY, the cancellation among diagrams is partial. The finite results obtained 
are related to the SUSY breaking scale Mgusy 

A naive expectation would be that Mgusy should be of the order of fiweak- This really has 
to be the case if SUSY is broken explicitly. However, if the supersymmetry breakdown is soft, 
then the splittings among supermultiplets are proportional to gMgusy, 9 being the coupling 
constant. This happens in supergravity: a soft breakdown may appear at Mgusy ~ 10^^ GeV, 
and mass splittings and M\y could still be of the order of 100 GeV. 

However, the 'small' scale, Mg^sy ~ 10^^ GeV, is still very large so that gravity cannot 
be neglected. Clearly, at such a scale all particles would have a gravitational interaction and 
the mass splitting would be of the order 



M2 

susy 



Planck 



10^ GeV, (1.12) 



i.e., comparable with the scale of electroweak symmetry breaking. The inverse Planck mass 
is basically the gravitational constant, k ~ Mp/^^^^. This leads naturally to local supersym- 
metry or supergravity. 

In addition to graviton (spin 2) there also appears gravitino, its superpartner (spin 3/2), 
with the mass of the order 

^3/2 ~ '"'^ , (1-13) 
^^J'Planck 

and the breakdown of supergravity induces the electroweak breakdown with M^y typically 
of the order of 7713/2. 
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2 Low-Energy = 1 Supergravity 



2.1 Soft Supersymmetry Breaking 

After the spontaneous breaking of a local = 1 supergravity one is left with a low-energy 
theory which is an explicitly broken global supersymmetry plus some 'soft-breaking' terms, 
the theory is the supersymmetrized minimally extended standard model. 

The full lagrangian of the SM contains superpartners as well as the standard particles. 
In the supersymmetric limit, the minimal model of supersymmetry has the gauge group 
SU{3) X SU{2) X U{1), with three generations of left chiral matter superfields 

Q(3,2,l/3), U^(3, 1,-4/3), D^(3, 1,2/3), L(l,2,-1) E^(l,l,2) (2.1) 

and two Higgs superfields 

H(l,2,l) and H'(l,2,-1). (2.2) 

The superpotential contains the most general gauge-invariant couplings which preserve mat- 
ter parity and contains only the fields required by supersymmetrization: 

W = h^UVrq^ + /i^H'D^Q,- + /i^H'E^^L,- + /xHH'. (2.3) 

The couplings hu, ho and He are 3x3 matrices in the generation space (i,j=l,2,3). 

The supersymmetric lagrangian is built in such a way that certain symmetries and con- 
servation laws are respected. For example, the baryon number B and the lepton number 
L are conserved. This automatically excludes terms such as (HL)^)^, (U'^D'^D'^)^)^, (LLE'^)^^, 
etc. 

The important new discrete symmetry introduced is it!-parity, defined as 

R={-f(^-^\-f, (2.4) 

where F is the fermion number. This symmetry is introduced to enforce baryon and lep- 
ton number conservation. With R-parity conserved, normal particles are R even, while all 
superparticles are R odd. As a consequence, the lightest superparticle (LSP) is stable. 

A global R invariance is achieved in the limit rfij — > and hu — ho — hE — 0, i.e., a 
chiral symmetry protects gauginos from getting mass in all orders in perturbation theory. 

The term /iHH' is introduced in order to avoid an unwanted massless axion. In the limit 
/X — > there exists an exact Peccei-Quinn symmetry, U{1)pq, which gets spontaneously 
broken when the Higgs fields develop vacuum expectation values. 

As we discussed earlier, the most convenient way of breaking supersymmetry softly is to 
couple it to a hidden sector of A^ = 1 supergravity. The situation is much simplified by the 
assumption of having a flat Kahlcr metric. There are a few types of soft SUSY breaking. 

(i) A cubic gauge-invariant polynomial in complex scalar fields (trilinear scalar couplings): 

[Cc/HU^Q + CdH'D^Q + e^H'E^L + a^^HH']^ + h.c, (2.5) 

where ^ are 3x3 fiavor matrices. This part of soft breaking contributes to squark and slepton 
mass matrices. 
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(ii) Gaugino Majorana terms: 



-Y,ma\aK + h.C. (2.6) 
^ a 

These terms obviously give masses to all gauginos, i.e., to wino, zino, and gluino. 

(iii) Mass terms of the scalar fields Zi of the chiral superfields: 

M^^z;zj + h.c. (2.7) 

These terms give masses to all scalar superparticles, i.e., squarks, sleptons, and higgs. 

(iv) The last term 

HBUW (2.8) 

enters the Higgs potential. 

In addition, if one assumes that the low-energy theory is a Grand Unified Theory, it 
would lead to a set of relations among the soft couplings and masses. Under the assumption 
that a grand unification appears at some scale Mx, a natural assumption is the universal 
gauge coupling. All soft-breaking terms are assumed to be parametrized at Mx by a universal 
gaugino mass mi/2, all scalar masses are given by mo at Mx, and there is a universal trilinear 
scalar coupling A. 

The above scalar and gaugino mass terms are allowed since gauge invariance is preserved. 
However, gauge invariance forbids mass terms for quarks, leptons, and gauge bosons - they 
will acquire their masses through spontaneous gauge symmetry breaking. Incidentally, only 
such particles have been seen until now. 

The only peculiar behavior as far as the origin of mass is concerned shows higgsino. While 
the higgs particle acquires mass from the usual weak spontaneous breakdown, the higgsino 
gets mass both from the /x-term in Csusy and from the weak SU (2) x C/(l) symmetry breaking. 

2.2 Flavor- Changing Neutral Currents 

The existence of two scalar doublets in the supersymmetric extension of the standard model 
does not lead to flavor-changing neutral currents since the underlying supersymmetry forbids 
it at the tree level. 

However, one-loop radiative corrections induced by the charged scalar would generally 
lead to the couplings between squarks, quarks, and gluinos containing flavor- changing parts. 
In addition, the hidden sector of supergravity, needed for the super-Higgs mechanism will 
bring new phase into the theory. If a typical SUSY breaking scale were of the order of 
10^° GeV, one might expect large renormahzation effects. 

Squark Mass Matrix 

After SU{2) x [/(I) symmetry is spontaneously broken, the potential invariant under 
global SUSY brings into play squark mass terms MqMq, Mq being the squark mass matrix 
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and with mass terms in £ typically of the form 



qlqL + {L-.R), (2.9) 

Ql,r being the squarks (L- and R- handed, respectively). In the minimal N=l supergravity 
extension of the SM, the soft-breaking terms contribute in such a way that squark masses 
are degenerate for different generations 

AqlQl + fJ'lqRqR, (2.10) 

and /i| ^ are of the order of the global SUSY breaking scale, /i ~ 10^ GeV, and usually taken 
to be equal to the gravitino mass 

Trilinear interactions between the Higgs field and two squark fields (cf. Eq. (|2.5| )) con- 
tribute with a strenght A and are proportional to fi. The parameter A may, in principle, be 
calculable from the hidden sector of supergravity. 

Left-right squark mixings, i.e., the mass terms of the type qLqR+ h.c. are induced through 
the electroweak spontaneous breakdown and are proportional to Mq. 

Finally, radiative corrections to the down-squark masses are induced by the charged Higgs 
and the higgsino loop and are proportional to Mj[M„, M„ being the up-quark mass matrix. 

The down-squark mass term becomes (taking into account quantum corrections) 

^2 ^ f /^il + MdMl + cM^Ml A*ms/2Md \ 

The parameter c is negative and is usually assumed to be 0{1). The parameter A is generally 
a complex parameter; it can be written as 

A = \A\e-^'^^. (2.12) 

The matrices M|, and cannot be diagonalized simultaneously by the same transformation. 
Let us introduce 2 unitary matrices, Ui and U2, such that U = U1U2 and 

f/i = ( 0^^ °^f/| ) ' ^^-^^^ 

where f/^, If^ diagonalize the down-quark matrix. Then the matrix U diagonalizes the down- 
squark matrix: 

WMlU = Ul ( + + ^^^^"^ l^l-^/^^^ ] U2. (2.14) 

Mu,d are diagonalized Mu,d matrices and 

V = Ul^Ut (2.15) 
is the usual Kobayashi-Maskawa matrix. 
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If one neglects the left- right squark mixing, and since ^ M|, a reasonable approxi- 
mation for U2 is 



which leads to 

^ [\A\7m/2MaV^ l^'l + Mj 
The current- and mass- eigenstates are connected via 



(2.17) 



d^U[ -Ifo) 1 (2.18) 
and the latter induce flavor-changing couplings. 
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3 Recent Approaches in Calculating Hadronic Transi- 
tion Amplitudes 

3.1 Chiral Perturbation Theory - Basics 

Chiral perturbation theory (ChPT) represents a viable alternative low energy theory of strong 
and electroweak interactions |^8|-[|83[. Its importance as an alternative approach to the usual 
formulation of the standard model is especially evident in hadronic processes where the lack 
of knowledge of QCD confinement is blurring our view of electroweak interactions. This 
is even more pronounced in kaon physics, since kaons and pions, being pseudo-Goldstone 
bosons, are even less reliably described in phenomenological quark models. 

This theory has recently become important in connection with some other alternative 

QCD simulation on latticed 



approaches, such as QCD hadronic duality sum rules ||56[1-[ 
and in the large- iVc approach of Bardeen et al. |^ . In all of these approaches, ChPT is used 



in some way or another: in duality sum rules, the tree level chiral realization of hadronic 
currents and/or operators is used in the parametrization of the hadronic side of the sum 
rule. In the large-Nc approach of ref. |7^, the chiral representation of currents is used to 
construct weak transition operators at zero momentum, which is needed in order to describe 
the operator evolution at large distances. Finally, the lattice calculation of weak matrix 
elements is (presently) performed for off-shell transitions, such as the K ^ tt transition. 
The relation with physical amplitudes {K vrvr) is achieved using the tree level ChPT 
relation between the above transitions. 

Chiral perturbation theory is based on our knowledge of the fundamental symmetries of 
the QCD lagrangian, e. g., the softly broken chiral SU{3)l x SU{3)ji symmetry. The formu- 
lation of ChPT is based USTl, |5^ on the folowing Ansatz: in any given order in perturbation 
theory, the most general lagrangian consistent with a given symmetry would result in the 
most general S matrix, consistent with incorporated symmetry and field theory conditions 
on anahticity, perturbative unitarity, etc. Formulated in this way, ChPT becomes a quantum 
field theory. 

The strong lagrangian at order (with minimal number of derivatives) is given by the 
nonlinear cr -model 

^trid^Ud'^U^) = ^gabi(l))d^rdy' (3.1) 



£(2) 



for a massless field, with the invariant metric 

gab{<P)=tr{dMdbU^). 



U is the unitary matrix field 



$ = :^A"0" is given by 



U{(f)) = expz-$. 



(3.2) 



(3.3) 



/ 71-° I rjg 

K- 



K+ 



(3.4) 



W ) 
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If the quark mass matrix Ai is different from zero, all terms in Cgtrong would pick up 
additional terms. Then, the strong lagrangian to the lowest order, reads 



42on, = ytr(a^f/9'^f/t) + V tr{MU + U^M), (3.5) 



with 

I = ^^1^ = = ^^1^ = Bo, (3.6) 

1 I — I (2) 

and V is proportional to the quark condensate, v = — ^ < 0|gg|0 >. The lowest order C strong 
at the tree level has basically two couplings, f.„ and which are expected to be calculable 
in QCD. 

The lagrangian Cgtrong of the nonlinear cr-model is nonrenormalizable and agrees with 
QCD only at the tree level (leading behavior) .0 In order to get the full strong lagrangian in 
ChPT one has to include all possible terms in the lagrangian and take account of all graphs 
in perturbation theory. The classical theory is equivalent to tree graphs of quantum field 
theory. However, a consistent quantum field theory requires the presence of loop graphs, 
since without loops, the theory violates unitarity. For example, the effective lagrangian 

£(2) 

describes the scattering of pseudoscalar mesons to order (tree graphs), the parameters 
being the quark masses, /tt and v. The unitarity demands that at order the T-matrix 
involves cuts with discontinuities, the contribution being determined up to a polynomial in 
the external momenta. In field theory, these discontinuities appear in the one-loop graphs 
from C^'^\ 

The inclusion of loops would generally lead to infinities. To get rid of them, the theory 
should allow for counterterms. Using the regularization that preserves chiral symmetry (e.g., 
dimensional regularization), one finds that one needs counterterms of order p^, since loop- 
graphs coming from the p^-lagrangian are of that order. Furthermore, one finds that all 
counterterms needed are contained in C^^\ This is intuitively clear, since we demand that 
regularization preserves chiral symmetry. Using the terms in as counterterms leads to 
finite results for all Green's functions to one-loop order. However, the renormalized graphs 
are unambiguous only up to a polynomial in the external momenta. 

In order to include the electromagnetic interaction, one has to couple the quarks to 



hermitian external fields, as proposed by Gasser and Leutwyler 
The lagrangian C^'^^ becomes 



o /■2 

£(2)= i^tr{D^UD^U^) + V tr{MU + U^M), 



(3.7) 



where is the covariant derivative 



D^ = d^U -ieA^[Q,U], (3.8) 

and is an external electromagnetic field. 

■^The choice of the effective lagrangian is not unique. At the tree level, both linear and nonlinear la- 
grangians lead to the proper value of the two constants and v. At the one loop level, however, they 
disagree. The linear (renormalizable) a model gives relations among the couplings which are in contradic- 
tion with experiment. 
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3.2 Chiral Realization of Hadronic Weak Interactions 



The extension of ChPT to weak processes is straightforward. However, the predictive power 
is often spoiled by a large number of unknown coupling constants. 

In weak AS = 1 kaon decays, CP violation, etc., one has only a few channels and this 
makes it difficult to go beyond the tree-level lagrangian. Nevertheless, the chiral realization of 
weak lagrangians||54|| has been found to be very useful as it 'supports' other approaches, such 
as QCD lattice calculations, QCD hadronic duality sum rules, and the large- A^c expansion. 

For example, to the order p^, the octet part of the AS = 1 weak lagrangian is given 
byH 

-^^aLi = 9s^s + hsQ, (3.9) 

where gs and hg are unknown coupling constants, not fixed by ChPT alone. The operators 
^8 and 9 have the following realization to the lowest order in derivatives and masses 

e = -^sicicg ^ tr{AUM + AiUM)^). (3.10) 
2v2 / 

The operator is the tadpole operator which in general contributes to the off-shell Green's 
functions, but does not contribute to the 5* matrix. For example, the operator 6 contributes 
to the amplitudes A{K tttt), A{K tt), and A{K —>■ 0). The first graph in Fig. l-I gives 
~ I ^sjma ^ which is , however, exactly canceled by the contribution of the second graph 
in Fig. 1-1. On the other hand, if one writes the amplitude A{K -kit) in terms of the 
amplitudes A{K —>■ it) and A{K —>■ 0) (the PCAC relation), one finds again that the tadpole 
contributions in the last two amplitudes cancel each other and A{K —>■ tttt) is tadpole-free. 



Recently, Shabalin|^ has suggested an interesting and intriguing mechanism that turns 
a quadratic GIM suppression of the d — s self-energy graph into a logarithmic one, show- 
ing a clear enhancement of the AI = 1/2 transition. The effect comes basically from the 
leading logarithmic one-gluon corrections to the bare graph. However, it has been shown by 
Guberina, Peccei, and Picek[^ that the full QCD correction in the leading logarithmic ap- 



proximation (LLA) reduces the tadpole contribution to a neghgible amount, which vanishes 
in the chiral limit. 

The physical amplitude A{K tttt) is therefore proportional to gs, which can be ex- 
tracted from experiment. Numerically, \gs\ ~ 5.1 . Next-to-leading corrections to iiT — ^ 2-7? 



and K ^ Stt have been calculated by Kambor et a/.||7^. The counterterms were fitted from 
the decay rates and slope parameters. The lowest-order value of gs reduces by 30%, whereas 
the analogous 27-plet constant remains practically unchanged. 



Effective Lagrangian for Radiative Kaon Decays 

(2) 

The electromagnetic field may be introduced into jCAg=i through a covariant derivative 
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^8 = ^^SiCiCs^s, (3.11) 

In addition, one has to add to it all possible terms of fourth order in derivatives and/or 



external fields allowed by chiral symmetry. Finally one finds [53 



em = --^F^'^iw, tr{QAL,L,) + tr{QL,kL,)} 

+ fIlk^^F''-F^, tr{AQUQU^) + h.c, (3.12) 

and Ln is a left-handed current with covariant derivatives. Only the terms relevant for 
radiative K decays with one pion in the final state are kept. Again, wi, W2, and W4 are 
undetermined couplings. To this order, one has to add the contributing part of the strong 
+ electromagnetic lagrangian 

C'^^rong+em = -teLgF^'tr {QD^UD.U^ + QD^U^D.U) + e^LioF^^F^^tr (QUQU^), (3.13) 

and the anomalous Wess-Zumino lagrangian 

^'Sz = ^^^e.^p.F^'^F'^^Tr' + ^) + O(0^ 0^), (3.14) 

linear in meson fields. This completes the tools necessary to calculate radiative K decays, 
which have been discussed in Chapter 4. 

3.3 Large- A/'c Expansion in ChPT and Beyond 



Large-Nc Expansion in the Weak Sector in ChPT 

The large- Ai'c limit can be taken consistently for all currents and/or operators appearing 
in the weak sector of ChPT. For example, the coupling in ( |3.9| ) is not determined in 
ChPT. The value of gg is known only in the large-A^^c limit, in which the weak lagrangian 
( |3.9| ) is reduced to the product of bare currents. With gs determined in the large- iVc hmit, 
the amplitude A{K nir) coincides with the large-A^^^ vacuum saturation approximation 
for the matrix elements of local 4-quark operators. 

The p^-level lagrangian introduces many unknown coupling constants. Typically, the 
amplitude receives contributions from one loop graphs calculated with the p^-lagrangian and 
from tree graphs stemming from the p^-lagrangian. The unknown couplings in p^-lagrangians 
are to be used as counterterms needed to renormalize possible divergences in loop integrals. 
The finite part of counterterms depends on the renormalization point /i; this dependence 
cancels the logarithmic fi dependence of the renormalized loop-graph contribution. If one uses 
dimensional regularization, one encounters only logarithmic divergences. Working instead 
with cutoff regularization, one also encounters quadratic divergences. The latter, however, 
disapppear after renormalization, as we shall discuss later. 
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The role of meson loops and/or counterterms is very important, as can be seen by studying 
rare kaon decays (cf. Chapter 4). The physical meaning of the counterterms (couplings) is the 
following: they are remnants of the original quarks and gluons after the latter are integrated 
out in the functional integral. Therefore, they contain both short- and long- distance physics, 
i.e., the effects of, e.g., vector mesons (long distance) and hard gluons (short distance). Of 
course, in some cases, such as the p'^ couplings Li,...,Lio in the strong+electromagnetic 
lagrangian, they are perfectly saturated by low-lying meson resonances. The situation is, 
however, quite different for weak couplings, e.g., the constants Wi in ( p.l3|) . Once determined. 



the renormalized couplings depend on the renormalization point /i; however, their natural 
order of magnitude is (47r)~^. 

The large-A^c expansion turns out to be a very useful tool in treating many problems 
in strong interactions. In the weak-interaction sector the situation is somewhat different. 
The idea of the laxge-N^. expansion ||75|| is based on the expectation that the true expansion 



parameter is not l/Nf. = 1/3, but rather something like l/AiiNc or even l/47iN^, as it 
happens in QED, where an expansion in the coupling constant e becomes an expansion in 
a = e'^/AiT. Unfortunately, in the weak-interaction sector the subleading 1/Nc corrections 
are often large, even huge in some case^. 



If one works with an explicit cutoff, the cutoff dependence comes in different ways||61 
(i) A possible contribution is absent because of chiral symmetry, (ii) quadratic cutoff 
dependence is, by power counting and chiral symmetry, of the form A^x tree-level C^'^\ (iii) 
logarithmic terms are of the form InA^x tree-level C^^\ It is then possible to absorb the 
cutoff dependence in the redefinition of coupling constants. It has been shown by Bijnens and 
Guberina[^ that the results obtained coincide with the ones obtained using dimensional 



regularization. 

As an example, the — mixing transition amplitude appears to be[]61 



< K'\C,,JK' >= (^ + ^)ml - ^5m^ + • • • (3.15) 

The leading term is of the order ~ i.e., it is the leading factorizable contribution. The 
new unknown coupling gi is of the order Nc- The cutoff dependence can be partially absorbed 
in the definition of the renormalized coupling f^'^"-^ and the rest enters the renormalized 
coupling gl^"'. The last subtraction is consistent with the large-Nc behavior, since gi and the 
quadratic divergence are of the same order in N^. The logarithmic divergences are handled 
in the same way. 

At the end, one can remove all quadratic and logarithmic divergences in the renormaliza- 
tion procedure. The logarithmic fi dependence is cancelled by the corresponding depen- 
dence of the counterterms. Both, the cutoff regularization and the dimensional regularization 



lead to the same result 61 



In the next subsection, we discuss the Bardeen-Buras-Gerard approach to the large-A^c 
expansion, and QCD hadronic duality sum rules. 



■^Tlic leading terms in the large-A^c expansion of the i^T-decay amphtudes lead to the ratio A{K'^ 
TT + TT^)/ A{K^ — > 7r"'"7r'^) of the order 1, which is by an order of magnitude smaller than the experimental 
result. 
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Bardeen- Bums- Gerard Approach versus ChPT 



The Bardeen-Buras-Gerard (BBG) approach ||71|| is based on the large-iVc expansion of 
QCD (with quarks and gluons), which imphes the existence of an equivalent dual- meson 
representation. The effective l/N^. expansion looks in the latter representation like a string 
theory. The leading- A^^c theory contains infinite trajectories of stable meson resonances. Since 
it is expected that the low-energy theory is only sensitive to low-lying states, one may use 
a proper truncation of the full theory. The proposal of Bardeen et al. is to use a nonlinear 
(T model, including loop effects, as a first approximation. The strong lagrangian (|3.5| ) in 
ChPT, discussed in preceding sections, now appears as the low-energy truncation of QCD. 

The correct infrared behavior of low-energy theory is guaranteed by low-mass and low- 
momentum loop contribution. The truncation is achieved by introducing a physical momen- 
tum cutoff in loop integrals. As expected, the main contribution comes from the quadratic 
cutoff A^. This is the crucial difference with respect to ChPT, where the cutoff depen- 
dence is absorbed in the renormalized quantities. Formally, the typical next-to-leading order 
amplitude in the BBG approach looks like 

2 

777 

c-\iiQCD)A ~ A^ + Mn + ■ ■ • , (3.16) 
which is to be compared with a typical amplitude in ChPT 



y4~w(/i') + 61n^ + ---. (3.17) 



For A ~ /i, the logarithmic terms are the same, and the finite counterterms in ( p.l7| ) are 



traded for the quadratic cutoff in (|3.16|) . However, this difference is only the formal one 
since the amplitude ( 3.171) is the full amplitude, whereas the BBG result ( p. 16 ) is the matrix 



element of the local operator, i.e., the short- distance behavior is factorized out in the form 
of the Wilson coefficient. This makes the comparison between the BBG approach and ChPT 
by no means trivial. In addition, the physical cutoff in the BBG approach plays the role of a 
scale that has to match the corresponding renormalization scale of QCD. In the chiral limit, 
for example, one has to match the quadratic A^ behavior with the logarithmic Iw^qcd one, 
which appears difficult ||6l]] . However, higher resonances are expected to improve matching. 



smoothing the quadratic behavior of the scale A and forcing the approximate logarithmic 
behavior[|. 

The second difference with respect to ChPT, i.e., that no counterterms were allowed in 



7T| means that, e.g., vector mesons and higher resonances have to be added separately. In 



ChPT, the latter are contained in counterterms. 

^ It is interesting to note[^ that one can exactly control the /i dependence if one calculates the leading 
1/Nc behavior of the penguin operator Qq = —^^^{s L<lR){<iRdL) ■ The leading term has the following chiral 
realization 



'i'^P^ = -mv^^tr [kd^Ud^V^). (3.18) 



Now, is proportional to {mq) ^, where niq is the running quark mass, whose QCD behavior is m? 



1 



rh^as^fMQQi))^''^^ ■ On the other hand, the leading 1/Nc behavior of the Wilson coefficient is given by 
c ~ as{fMQ(j]j)^^^^ ■ Obviously, doing a systematic l/Nc expansion, one achieves an exact cancellation of the 
fi dependence. 
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As we mentioned, one would naively guess that the quadratic cutoff would have to mimic 
counterterms if both theories were equivalent. This, however, is not true for the following 
reasons: (i) contrary to the cutoff-terms, the counterterms contain long-distance effects 



as well as short-distance ones. In the approach of [^, the latter are factorized in the 
form of Wilson coefficients; (ii) the counterterms in ChPT also contain the effects of higher 
resonances (vector mesons, etc.). As given in ||71|, they have to be added separately. 

In addition, ChPT would implicitly require the importance of counterterms in order 
to keep the whole correction in weak amplitudes moderate. Otherwise, the perturbative 
expansion would break since in some cases loop corrections contain huge logarithms ||76[ . 



This implies the importance of missing corrections due to vector mesons, etc., in the BBG 
approach. 

From the above points it appears difficult to achieve one-to-one correspondence between 
ChPT and the large-A^^c approach of Bardeen et al. However, we would like to point out 
that in spite of our inability to achieve it, this genuine approach of Bardeen et al. leads to 
reasonable results. For example, the result for the B parameter, B ~ 0.7 is rather stable 
as well as the result for the — > tt+tt" amplitude. The latter is actually close to the 
experimental value, i.e., the Al = 1/2 rule seems to be explained in the BBG approach. 
The Al = 3/2 amplitude, however, is sensitive to the choice of the cutoff which induces a 
large uncertainty. 

Recently, an exciting project has been started|7^ in an attempt to derive the low-energy 



chiral realization of the AS = 1,2 operators by integrating out the quark fields in a gluonic 
background. The spontaneous breaking of chiral symmetry is triggered by a source term, 
which is added to the QCD lagrangian. The bosonization of the QCD lagrangians generates 
the effective lagrangian of ChPT but with explicit values for different couplings! Further 
work in this direction is in progress. 

3.4 QCD Hadronic Duality Sum Rules 

In the standard model, the effective hadronic weak lagrangian is usually given in terms of 
local 4-quark operators, e.g., the AS* = 1 lagrangian is given as 

'^e// = Ec.(/^)(^«, (3.19) 

i 

where Cj(/i) are the Wilson coefficients and Oi are local operators. As we have discussed in 
the preceding sections, there is a complementary representation in the form of effective chiral 
lagrangian ( p.9| ). Both pictures are equivalent. From the technical point of view these are 
associated with definite problems: using the lagrangian ( p.l9| ), one cannot reliably calculate 
the matrix elements of the composite local operators and the use of the complementary in 
ChPT gives rise to a large number of undetermined couplings. The duality proposed in the 
QCD hadronic duality sum-rule approach spells out the consistency of both representa- 



tions; indeed, there is a 'window' where both pictures reliably describe the same quantity. 
This is achieved by writing down a system of finite-energy sum rules (FESR's) which relate 
the integrals of the hadronic spectral functions to the corresponding integrals calculated in 
QCD. 
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The starting point is the two-point function 

^(g2) = i J d^xe'''-^ < O|T(O(x)0t(O))|O >, (3.20) 
whose imaginary part enters the sum rules 

K„(So) = / dte -Im ij{t)hadronic = W / rft T ^ | < 0|£|r > 1^ 
Jo 7T Jo p 

Sq^^ as{so) ,as, - ,m2(so) 



{a + a' — (so) + h- 



(n + 5)(167r2)3 ^^(/i^) vr Sq 

m, < OlgglO > < Ol^s^F^^F'^^^jO > ^ ^ ^ 

+c— + rf ^^^^ + ...}, (3.21) 



where g' = c ^{jJ^g. As can be seen from the Lh.s. of ( |3.21| ), the spectral function ^Imipit) 



derived from the correlator (p.20|) describes, how the operator O couples the vacuum to 



physical states. In its ChPT version, it enables one to calculate the threshold behavior of 
the intermediate states (F = Kti, Ktctc , etc.). The unknown constant g' is factorized out 
as an overall normalization. The upper limit of integration, sq, is the onset of the QCD 
continuum. It should be high enough, so that the QCD expansion on the r.h.s. of ( p.21| ) has 
sense. The window we are talking about is roughly of the range of Sq where the two pictures 
are dual. Clearly, this window cannot be of large range; for too low values of Sq, the QCD 
expansion breaks, and for too high values, ChPT cannot reliably parametrize the hadronic 
spectral function. 

The parameters a and a' on the r.h.s. of (|3.21|) are calculable in perturbative QCD; a 



gives the asymptotic behavior in the chiral limit, and a' is a coefficient of finite corrections. 
Mass corrections are taken into account by the second term, and the third and fourth terms 
are nonperturbative corrections coming from quark and gluon condensates, respectively. 
Leading logarithmic corrections are summed up by making use of the renormalization group 
equation, and they are factorized out as a ratio of coupling constants to the power of 7, 
the latter being proportional to the anomalous dimension of the local operator. For the 
multiplicatively renormalizable operators, c^(/i^) is of the form 

c\^') = (^^)-^ (3-22) 

i.e., the /i dependence of the Wilson coefficient exactly cancels the fi dependence of the 
matrix element. 

In order to make sense, the power corrections in ( p.21|) should not exceed 20 — 30% 







otherwise, higher mass and/or condensate corrections become important. This basically 
determines the allowed range of Sq from the point of view of perturbative QCD expansion. 
It turns out that in the — mixing [B parameter) [^] , and in the tt^tt^ decay 

(A/ = 3/2 transition) ||57||, this range is rather high, so ~ 7 — 11 GeV^. This implies that 
one has to correct the pure ChPT behavior of hadronic spectral functions for the formation 
of resonances. The proposal in is to modulate the final-state interaction through the 



Breit-Wigner form and normalize it in such a way that for t = it reduces to the chiral limit 
value. 
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The onset of scaling is obtained by taking the ratio of the two sum rules 



^n^O ~ CO . s , (3.23) 

so that the unknown constant g' drops out. The ratio r is usually normalized in such a 
way that that it equals 1 for the asymptotic case (no QCD). Then one looks for the duality 
region, i.e., for the region where the function r„(so), calculated in ChPT, agrees with the 
same function calculated in QCD. Once the duality range of Sq is fixed, any of the sum rules 

leads to the value of g, for any sq in the duality range. With g determined, one can easily 
get the decay amplitude using the ChPT lagrangian. 

The B parameter plays the role of the constant g' . It has been calculated by Pich and 



de Rafael p6|, with the result 

\B\ = (0.33 ± 0.09)[a,(/i2)]2/^ (3.24) 

Compared with typical results obtained in the alternative approach (large- Nc expansion, 
lattice QCD) the result appears smaller roughly by a factor of 2, and, obviously, the 

vacuum saturation larger by a factor of 3. 



The same method has been applied by Guberina, Pich and de Rafael ||57[| to the K 
7r+7r° decay and the result ( |3.24| ) agrees with experiment within a few percent. This should 
be considered as a successful test of the duality approach. It should be noted that especially 
this amplitude has almost always been reproduced in different approaches within a factor of 
2, but better accuracy is difficult to achieve. 



Recently, Prades et a/.|^ have improved the calculation of -B by a careful analysis of the 



hadronic parametrization. The result is somewhat higher, \B\ = (0.39 ± 0.10)[as(yU^)]^/^ 

Unfortunately, the success shown by the above results does not pertain to the A J =1/2 
amplitude. The calculation by Pich et al. failed to reproduce this amplitude by an order 



of magnitude. It was quite difficult to understand this failure. Recently, definite progress has 
been made|S^ which points out the reason for the failure. Namely, in the calculations of the B 



parameter and the A/ = 3/2 amplitude, finite as corrections were taken into account. They 
were found to be moderate, and because of technical complexity, they were not calculated 
for the A/ = 1/2 amplitude, i.e., it was assumed that they were also moderate in this case. 



However, Pich has recently shown |]62| that this assumption is premature. The perturbative 
finite as corrections to the A/ = 1/2 amplitude are so huge that the whole perturbative 
expansion breaks. This clearly shows that the problem is highly nonperturbative, and 
cannot be handled in a perturbative way. 

The calculation performed in includes two assumptions that significantly reduce the 



complexity of calculation: (i) The operators 0± in the bare weak lagrangian are handled 
without mixing, i.e., as they are multiplicatively renormalizable, (ii) the penguin operator 
Qq is taken in the large-A^^c limit in order to be multiplicatively renormalizable. Then, the 
calculation shows that for finite as corrections are moderate, and for and Qq, they 
appear with the coefficients 47/5 and 423/20, respectively. With the usual /i^ = t rescaling 
which eliminates all logarithms in the spectral function, the corrections exceed 100%, even 
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at t = 10 GeV^, showing a clear breakdown of the perturbation expansion^. 

We would like to point out that the discussed breakdown of the perturbative QCD ex- 
pansion for AS* = 1 transitions may have serious consequences. Clearly, it is not possible to 
use the QCD duality approach in calculating these transitions. However, both the laige-Nc 
approach of Bardeen et al. and the lattice calculation also start from the Wilson expansion 
with the fi scale set at hadronic level. The huge radiative corrections both in the Wilson 
coefficients and in the matrix elements cast serious doubt on the validity of the starting point 
in the above calculations. 



3.5 Hadronic Matrix Elements on the Lattice 

Lattice QCD offers a method for calculating hadronic matrix elements from first principles. 
The technical difficulties are however enormous. Besides the usual problems caused by 
insufficient lattice size, lattice spacing, and computer time, the matching of the lattice and 
continuum operators is complicated by a conflict between the lattice regularization and the 
chiral properties of the theory [^. Since the chiral symmetry is broken on the lattice with 
Wilson fermions and the flavor content of the standard model is broken on the lattice with 
staggered (Kogut-Susskind) fermions, one must take a linear combination of lattice operators 
in order to form the desired continuum operators. This 'mixing problem' is extensively 



discussed in ref. []64 |. 



The lattice calculation of the matrix element of a local operator between meson states is 
based on 'measuring' the corresponding three-point correlation function on the lattice: 

G{x, 0, y) =< P,{x)O{0)PM > . (3.26) 

P5 are diquark operators with the same flavor content as the corresponding mesons, and O 
is a 4-quark operator, e. g., appearing in the effective lagrangian (|3.19|) . Let us illustrate 
the method on a typical matrix element, for example < 7T^\sruurd\K^ >. By replacing 7r+ 
and by their diquark operators we arrive at the correlation function 

0, y) =< rf(x)75M(x)s(0)rM(0)M(0)rd(0)M(2/)75s(y) > • (3.27) 

The contraction of the quark flelds in ( pj.27D in all possible ways yields the 'eight' and the 
'eye' diagrams depicted in Fig. 2. In K vrvr, the eye diagrams are pure AI = 1/2, 
whereas the eight diagrams may be either AI = 1/2 or 3/2. In the AS = 2 transition, only 
the eight diagram contributes. The propagators in Fig. 2 are quark propagators averaged 
over gauge conflgurations on the lattice. The conflgurations which include internal fermion 
loops are absent if the so-called 'quenched' approximation is used. Technically, the fermion 
determinant in the averages is set equal to one. The approximation is justifled due to al/Nc 
suppression of the femion loops. In the eye diagram, however, the w-quark loop is not of 



^ It is interesting to note that this behavior persists in in higher orders. In the large- A^c hmit, the 0{ag) 
corrections can be easily computed, leading to the spectral function of the penguin operator Qq: 

mpeng ~ [«. (i)] ^^^^ { 1 + 24.28^ + 470.72(^)2 + ...}. (3.25) 
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higher order in 1/Nc and therefore the eye cannot be ehminated on the basis of the quenched 
approximation. Moreover, its contribution is beheved to be responsible for a large part of 
the A J = 1/2 enhancement. 

As regards the K ^ Tin decays, there are basically two methods of calculation, depending 
on the prescription used for lattice fermions. In the first method, the K — >■ Tin amplitude 
is calculated directly, with all particles at rest and all quarks degenerate. All particles are 
on-shell, so the amplitude is well defined but the weak Hamiltonian must insert momentum. 
Lowest-order ChPT is then used to relate the amplitude to the physical one. The method 
involves a relatively small chiral extrapolation and Wilson fermions can be applied. This 
method has been used by the Bernard and Sony group p3, [T3, B6| and the European Lattice 



Collaboration (ELC) H, [T|, |67 



The second method is to measure the K ir and K vac matrix elements and use 
ChPT to relate these to the physical K ^ -ktx amplitude [Q]. In ChPT, the physical ampli- 



tude is proportional to the matrix element < 'KlOsublK >, where the subtracted operator is 
defined as 

Osub = 0- pirns + md)sd. (3.28) 
The coefficient p is determined from the ChPT relation 

< 0\O\K >= p{ms - ma) < 0\s-f5d\K > . (3.29) 

This method is appropriate for staggered fermions owing to their good chiral properties, 
which allow the use of equation ( |3.29| ). In addition to calculating the eight and the eye 
diagrams, one has to calculate the matrix elements of the two-quark operator subtraction 
and the K vac diagram needed to fix the coefficient of the subtraction. This method has 
been mainly used by the staggered group [^, . 



If one computes K ^ tt or K ^ Tin on the lattice, one has to face up to dealing with the 
TTTT interactions in the final state (FSI). The FSI may cause a significant enhancement of the 
AI = 1/2 amplitudes, and suppression of the AI = 3/2 amplitudes |6^. In addition, there 
is a final state phase to be determined in order to sum both amplitudes to get the physical 
K —>■ TT+TT^ rates. 

In the last few years calculations of the hadronic matrix elements on the lattice have 
shown substantial improvement. The systematic errors and finite size effects are now under 
better control. Let us quote some of the recent results most of which, however, should still 
be considered as qualitative. 

Relatively reliable data are obtained for the KK amplitude. The ELC collaboration 
quotes two values for the B parameter at a~^ = 1.34 GeV, depending on what type of fit 
is used: 0.91 ± 0.11 (linear fit) and 0.64 ± 0.11 (quadratic fit) The staggered group 



observed a deviation from the asymptotic scaling in the region of (3 = G/g"^ from 6 to 6.4. 
They find that if the lattice spacing is reduced by a factor of 2 (keeping the physical volume 
fixed), B drops from 0.69 ± 0.02 to 0.54 ± 0.05, whereas the scaling requires a change of 
about a few percent. The extrapolation to a = with the anomalous dimension included 
would give Bg~^^^ 0.45 1|6^. Bernard and Sony observed a similar behavior. Going from 
the lattice size 16^ x 40 to 24^ x 40 at /5 = 6, they found that B dropped from 0.86 ± 0.24 to 



0.69 ± 0.12 A probable explanation for this scaling violation is that one still sees large 
0{a) effects. 
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As regards the K ^ tt and K ^ tttt amplitudes the data are less reliable. The ELC 
group used two methods of calculation. The first one was to translate their KK data into 
the Ktt amplitude and the second was a direct — > vrvr calculation. The results for the 
AJ = 3/2 amplitude obtained using the two methods are = (7.0 ± 0.8)10~^mi^ and 
(8.6 ± 0.8)10~^mK, respectively, to be compared with the experimental value 3.7 ■ IO'^itik- 
The ratio R obtained using the K ^ n method is 12 ± 5 with 75 configurations, but 
unfortunately 12 ± 12 with 110 configurations. The direct method yielded a worse result: 
i? = 35±30 [0. 

Sharpe suggested that the AI = 1/2 rule might be an accumulation of factors of 1.5 - 
2 due to different mechanisms ||69[. Taking into account an enhancement of R arising from 
the FSI, the lattice calculation should aim at about 62 % of the experimental value, i. e., 
Reye+eight ~ 14. The staggcrcd group finds that the eight diagram alone gives Reight = 
3.6 ± 0.02, but according to their last data, the eye contribution is unexpectedly consistent 
with zero ||7^. 
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4 Rare Kaon Decays 



4.1 Kl and K^^s 77 Decays 

As far as weak interactions are concerned, the Kl — * proceeds partly via a box diagram 
(diagram (a), Fig. 3), and partly via an sdZ vertex (black box in diagram (b), Fig. 3). The 
first diagram is similar to the process — > 77, which proceeds via diagram (b), Fig. 3, 
where q denotes up-type quarks, u, c, and t. In fact, a major contribution to the decay rate 
for Kl — >■ fJ'fi comes from the higher-order electromagnetic process (digram (d). Fig. 3), i.e., 
via a two-photon intermediate state. This diagram actually dominates the imaginary part 
of the amplitude A{Kl — >■ /x/i), relating it to the process K^ —>■ 77: 

-r^absorptive -n / a -i \ 

^KL-,t^il ~ii^L-77 (4-1) 



We would like to stress that eq.( [4.1| ) would, in principle, provide us with a lower limit for 
the partial width 

-f^absorptive /a qx 

with being the 'unitary bound'. 

Although, naively, the processes K^ — > /U/i (diagram (a)) and K^ — > 77 (diagram (b)) 
would both have amplitudes of comparable strength, A ~ jM^^ ~ GpO--, experimentally, 
they differ by a factor of 10~^, the reason being, of course, the GIM suppression. The 
respective branching ratios are given bylB 



BR{Kl fifi) = (7.3 ± 0.4 X 10^^ (4.3) 
BR{Kl 77) = (5.70 ± 0.27) x 10"^ (4.4) 

Closely related to the Kl /i/i decay is the process K~^ — > n'^ui'. It receives contribu- 
tions both from diagram (c) and diagram (e) in Fig. 3. Especially, the unitary bound derived 
for Kl would set constraints on the tt^uu decay, provided the Kl fifi decay 

rate were dominated by the absorptive part of Kl ^ 77 ^ /i/i and the dispersive part were 



negligible [80 



The effective lagrangian for the ds — > fifi and ds — > processes in zeroth order in strong 



interactions is |79| 



(^2 /M'2 

C = ^^snMC^Lrf^L (4.5) 

^ i 

The coefficients C and Di (i = generation index) are related, via the unitarity relation 
'^js^jd ~ ^1 Inami-Lim coefficients C and -Di|7^ which are functions of the mass 

ratios Xj = rn^^jMy^ and yi = mfjMyy. 

As we discussed earlier, there are potentially large long-distance effects in the Kl — ^ fifi 

decay. 

It is usually assumed that the dispersive part of the amplitude for Kl — > 77 — > is 
small compared with the absorptive one|5^. If this were true, the limit on the short-distance 



part of Kl — * /U/i would be obtained as 

A{Kl fifi)short-dist. < A{Kl fip,)exp - A{Kl ^ 77 ^ ^ijl) absorp.- (4.6) 
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However, if the dispersive part of A{Kl — 77 ^ /i/i) were large, and with opposite sign to 
A{Kl — > fJ'fi)short-dist., the above constraint would no longer be valid. 



The decay jjfl has been studied in ChPTpiI]. The main contribution comes from 

the diagram (a). Fig. 3, i.e., the decay proceeds via the 77 intermediate state. The P77 
vertex is described by the Wess-Zumino term. The corresponding integral is logarithmically 
divergent and the theory requires counterterms in order to renormalize the divergences. The 
finite part of the counterterms unfortunately cannot be determined, and one relies only on 
the logarithmic terms. This is certainly not quite reliable, although one does not expect that 
this arbitrariness essentially changes the conclusions. 

If one performs renormalization in the MS'-scheme, the two-photon intermediate-state 
contribution is||8T[] 

r^.^., = r^.^..^(^)VP, (4.7) 

where Re A contains the logarithmic term, which is, of course, dependent^ : 

«e A ^3Hf) - 7 + i{iln^(l^) - H- 2Lm[^)}, (4.8) 

and P = (1 — Amj^/mj^Y^'^. On the other hand, Im A has no /i dependence: 

/m^ = -^ln(i±|). (4.9) 
The amplitude A has been obtained by using a once-subtracted dispersion relation with 



Im A being known from the calculation performed in ref . [p2 
Choosing /x^ ^ one gets 

^^^^^ = 3.6 X 10-^ (4.10) 

Theoretical result is about twice as large as the experimental one. This, however, should be 
considered as reasonable agreement, since: 

(i) the calculation is not complete; the finite part of the counterterms has not been 
determined, and 

(ii) the /i dependence of Re A is pronounced. In addition to the rather strong /i de- 
pendence (a factor of 4 in the range 0.2 < < 1 GeV ), one finds that, starting from 
/i ~ 0.3 GeV , a dispersive part dominates over an absorptive one, the former being a factor 
of 2 larger than the latter, for /i^ = m\. 

Having in mind points (i) and (ii), it would be premature to claim the real dominance of 
the dispersive part. However, contrary to the usual assumptions, this shows that it may be 
significant. 

Before we compare this calculation with the calculations in different approaches, we 
should remember the following. In principle, if the counterterms were determined, this would 
be a complete calculation of F^-^^^^; i.e., it would also contain short-distance contributions. 

^ If one were able to determine the finite part of the counterterms, the /i dependence would disappear, 
i.e., the counterterms are also /i dependent. 
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In addition, all effects of heavy particles appearing in the theories which go beyond the SM 
would also be included in the counterterms. 



Ks — > 77 process in ChPT 



This process is described in ChPT by the lagrangians (|3.11| )- (|3TT^ ). In the calculation 
performed in ref. only the A J = 1/2 part in £as=i has been kept. Then, the process 
proceeds via the diagrams shown in Fig. 4. It is the simplest example of the predictive 
power of chiral perturbation theory because it can arise only from loop graphs involving 
but does not receive any contribution from the direct couplings wi, W2, w^. These 
couplings play the role of counterterms, i.e., the loop divergences have to be absorbed in 
these couplings. Since couplings are absent, that means that divergences must cancel and 
loops have to be finite. The graphs displayed are expected to have both quadratic and 
logarithmic divergences. The former are cancelled by an SU{?>) invariant counterterm, and 
the latter also disappear. This surprising result is a consequence of the requirement that the 
amplitude vanishes in the m\; = limit. The final result is given as|]8l 
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(4.11) 



where 



F{z) = {1--[tx^ - hi^Q{z) - 2m\nQ{z)]}, 

l + (l-4/z)V2- 



(4.12) 



The amplitude entering ( |4.11| ) is dominated by the imaginary part, which is also given in 
ref. [^. Numerically, eq. ([4.11|) gives 



1.4 X 10^2° GeV. 



From ( ^.13| ) one infers that 

BR{Ks ^ 77)theory = 2.0 X 10-^ 
The measured branching ratio is|TI5] 

BR{Ks ^ 77)exp = (2.4 ± 1.2) x 10-^ 



(4.13) 



(4.14) 



(4.15) 



in excellent agreement with the theoretical prediction ([4.14|) . This is a nice example of how 
the full ChPT completely accounts for the process in question 
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Kl 77 in ChPT 



This decay basically proceeds via the vr^ , r], and rj' poles 



^ AiP.^-n)<P.\n\K,>^ (4.16) 



The rj' contributes only through an octet weak transition. The vertices P77 are again 
described by the Wess-Zumino term. The amplitude (|4.16|) has been carefully calculated in 
and by Goity [Bl]. There are two sources of uncertainties that enter that calculation. 



The first one is the rj — rj' mixing angle 9, whose previous value was about 10°. However, 
recent measurements of rj production in 77 collisions[|^ give a higher rj ^ 'y'j rate, and 
consequently a larger mixing angle The rj and rj' decays are well described by taking the 
values [p3l 

= 24° ±2° and /^//;) = 1.04 ± 0.05. (4.17) 

In spite of that, the result for A{Kl 77) is very sensitive to the parameters 6 =< K^lrjg > 
and K, =< K^lrji >. Plotting the experimental BR{Kl 77) in terms of 5 and k shows 
that for agreement with experiment, the parameter 5 should deviate from its chiral limit 
value0. 



4.2 Radiative Kaon Decays in ChPT 

We have already discussed the decay Ks — > 77, which is a nice example of the self-consistency 
of ChPT. This decay belongs to class (i) which includes decays where the contributions from 
the counterterm lagrangian are forbidden by some symmetry. Class (ii) includes decays where 
loops are finite as in class (i), but there is an additional contribution of the counterterms 
that is now scale- independent. Finally, class (iii) includes decays where the loop amplitudes 
diverge. In this case, ChPT should allow for the counterterms. 

Decays K —>■ 7777 

The amplitude n'^'-ff can be completely calculated in ChPT. It is uniquely determined 

since, as in the process Ks — > 77, the counterterms do not contribute, and the amplitude is 
given in terms of the loop amplitude. The spectrum is given by[Q 

dy " ~(4^0^ )\yy-^ )F{yz) + {l-y + z )F{y)\ 

y = m^^^/ml [0 < y < (1 - z)-^ = 0.52] 

A(a, b, c) = a^ + 1)^ + - 2{ab + bc + ca), (4.18) 

® A large mixing angle, 9 « 20°, is also predicted in the large- A'^c limitp3] 

^ The good agreement with experiment, obtained for 9 ~ 13° in ref. |pl[, is, of course, spoiled with the 
present value of 9. The same is true for the calculation of ref. [p5[. 
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where F{z) is the function given in ( [4.12| ). The spectrum has a characteristic shape, quite 



different from the phase space. The integrated rate (|4.18| ) gives 



BR{Kl 7r°77) = 6.8 x 10"^. (4.19) 



This branching ratio has recently been measured. The reported result |T9[| is: (2.0±0.5)xl0 ^. 
The parameter free prediction is the ratio 

^^^^^ = B.9 X 10- (4.20) 

T[Ks 77) 

A similar process, —>■ 7r+77, belongs to class (ii): contributions come from a finite- 
loop amplitude, from scale-independent counterterms|3j and, in addition, from the anomaly. 
The lower bound on the branching ratio is predicted to be 

BR{K+ 7r+77) > 4 x 10~^ (4.21) 



and is expected to be measured soon, since the present upper limit is [p^ 

BR{K+ 7r+77) < 10"^ (4.22) 

The decays K — > 7177 are also interesting because CP violation appears owing to the 
interference between the absorptive amplitude and the counterterms (complex numbers). 
The charge asymmetry is given as 

T{K^ 7r+77) - T{K- 71^77) = Im c 1.5 x 10"^^ GeV, (4.23) 

which crucially depends on the value of c. There is also a charge asymmetry for the K 
7ie~^e~ decays; it is given as 

T{K+ TT+e+e") - T{K- 7r-e+e") = Im w+ 1.6 x 10"^^ GeV, (4.24) 

where is the following combination of counterterms 

w+ = -^^{wi + 2w2-l2Lg). (4.25) 

The source of CP violation in the standard model is presumably a phase in gg. Turning on 
the electromagnetic interaction induces the electromagnetic penguin as a new source of CP 
violation. The latter contributes only to the counterterm Wi in the large- iY^ approximation. 
Since the counterterms Lj are real, the following relation holds: 

327r2 

Im c = 2 Im = Im Wi. (4.26) 



This relation predicts the ratio of charge asymmetry ( |4.23|) versus charge asymmetry ( |4.24| ) 
to be ~ 200. 



The counterterms enter the amplitude in the form c — 327r^[4(L9 + Liq) — ^{wi + 2w2 + 2^4)], which 



is estimated to be 0(1) [p3|. 
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Decays Kl vr'^e+e and rj 7r°77 



The one-photon exchange contribution to the amphtude for vr^e^e" is purely CP 

violating. The intrinsic CP violation is again related to Im wi and is comparable with the 
'normal' CP violation. However, what makes this amplitude peculiar is the rate: the one- 
photon exchange leads to BR ~ 10"^^ — 10~^^. Higher-order contributions in ChPT|53] are 
by two orders of magnitude smaller than the one-photon exchange contribution]^. 

The electromagnetic suppression of the CP-conserving amplitude makes it probable that 
this decay is dominated by the CP-violating contributions coming from the small CP-even 
Kl component of the and through the direct CP- violating contribution in K2 — > 7r^e~^e~ . 
The latter is of the same order or even larger then the indirect one. 

There is another candidate, the rj — > 7r°77 decay, which causes difficulties. It receives 
contributions only from the finite loops, and the unique prediction is 

T{r] ^ TT^-f-f) one-loop = 0.35 X IQ-^ eV. (4.27) 



The experimental value is larger by two orders of magnitude [[19 



r(r^ ^ A7)exp = (0.85 ± 0.19) eV. (4.28) 

Higher-order contributions are estimated, and the result is still far from the experimental 
value. This decay really looks like an unexpected failure of ChPT; however, before claiming 
it, one should rather wait for confirmation of ([4.28|) in independent measurements. 



As expected, correct introduction of vector mesons as nonlinear realizations of chiral SU{3)v leads 
basically to the same results as those obtained using counterterms. The blind use of vector-meson dominance 
produces, e.g., contributions which are not permitted by chiral symmetry[^. 
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5 Rare Decays n^uu and 7r+77 

5.1 Supergravity Effects in Rare Processes 

Given the fact that the hmits on superparticle masses are presently rather high, roughly of 
the order of M\y or higher, it is clear that supersymmetric effects are expected to be very 
tiny. The strongest constraints seem to come from the /i — > 67 process in the lepton sector 
and the — Ks system in the quark sector. The constraints imposed by these two processes 
are in general respected in — > /i/x, — > fie, jj, eee, nN — >• eA^, etcPO, pT 



Muon anomalous magnetic moment 



No useful constraints are expected from the g — 2 factor The agreement between 
theory and experiment for the anomalous magnetic moment of the muon is better than 
2 X 10^®. In the supersymmetrized self-energy graph of the muon, the photino and the 
smuon enter the graph. With fi^ and fiR degenerate, the contribution is proportional to the 
square of the muon mass. Then, the effective vertex takes the form 



2m^ 



F{q )ua^yq''u, (5.1) 



where q is the momentum of the photon. The g — 2 gets a contribution through -F(O) which 
is proportional to the muon mass. Besides, there are graphs with wino and zino. Fig. 5. A 
massless photino graph has been calculated by Fayet [p^] with the result a = l/2{g — 2) ^ 
10~^, provided the smuon masses are larger than 15 GeV. Adding wino and zino graphs. 



Ellis, Hagelin and Nanopoulos|p3||, and Barbieri and Maiani[^ have calculated the full 
contribution in global SUSY, with essentially the same result. The extension to a more 
realistic supergravity model has been made by Romao et al.[^. Contributions of neutral 
gauginos have terms quadratic and linear in m^. The latter are potentially large, but they are 
suppressed owing to a kind of GIM mechanism[^]. At the end, the dominant contribution 
comes from the photino graph, which dominates by an order of magnitude. For a given 
photino mass, a is a decreasing function of gravitino mass. For 7/13/2 < 15 GeV, the present 
agreement between theory and experiment is preserved. This means that supergravity effects 
on g — 2 are really small. 



H e'j decay 



In the standard model with massless neutrinos the decay ^ 67 is strictly forbidden 
owing to lepton flavor conservation. This also remains true in the supergravity model. 
However, allowing for massive neutrinos compatible with present experimental upper limits 
on their masses, the branching ratio for this decay is of the order of ~ 10"^^ in the standard 
model. This is far below the experimental upper limit, which is 7.2 x 10"^^. 
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In the supergravity there is an additional diagram for the /x — 67 decay, diagram (a) in 
Fig. 5. 

Keeping the same values of neutrino masses and mixing angles as in the SM, and varying 
gravitino and photino masses in the range 20 — 250 GeV, one finds that for a given 7713/2 
the branching ratio increases with m^j,, the reason being that the increase of the photino 
mass results in the decrease of the wino mass. Since the branching ratio contains a term 
proportional to (M^y/M^)^, it produces an enhancement for light winos. This enhancement, 



however, weakens with increasing gravitino mass|^5|. Even for a very light gravitino, 7713/2 = 
20 GeV, the SUSY branching ratio is smaller than 10^^^. 

In view of the present experimental limits on SUSY particles, these effects are also very 
tiny. 



Kl — Kg System 



In the standard model, the Kl — Ks system is described by box graphs with u, c, t quarks 
entering the loops. For degenerate quarks, these contributions cancel. The same remains 
true for supersymmetrized graphs since the partners of quarks and W bosons, squarks and 
gauginos, have the same quantum numbers. 

This cannot lead to the limit on squark masses, but gives limits on mass differences 
between the families. The bounds obtained are even strenghtened if one includes graphs with 
strong interacting gluinos. The condition for this is that superscalars have an off-diagonal 
mass matrix in the basis where the quark mass matrix is diagonal. Thus, the bounds on 
masses of d, s, b superscalars are obtained. These bounds, coming from gluino exchange, 
are generally by an order of magnitude stronger than the corresponding ones coming from 
simple supersymmetrized graphs. 

The usefulness of bounds is often spoiled by uncertainties in the QCD calculation of 
matrix elements of the effective operators once the heavy (super)fields are integrated out. 
We discuss these problems in Chapter 6. 



5.2 The decay tt^uu in the standard model 

At zeroth order in strong interactions the decay -k^vv proceeds through diagrams 

in Fig. 6. The first group of graphs, group (a), are box graphs, which, calculated in the 
i?g-gauge, contain a gauge- dependent part 7(0- The same function, 7(0) also appears in 
group (b) of diagrams, but with an opposite sign. Therefore, the whole set of graphs in Fig. 
6 is gauge independent, as it should be. The whole set of graphs in Fig. 6 leads[^, ^ ^ 
to the result ( |4.5| ) for the effective lagrangian. 

The branching ratio BR{K^ n^uu) can be related to the BR{K~^ TT^e'^v) by using 
an isospin relation [[7^] 

BR{K+ n^v.v,) = 77 7/2' BRiK+ n'e^u,). (5.2) 
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Using the experimental valuesQ BR{K+ n^e+Ue) = (4.82±0.06)%,Mh/ = (80.22±0.26) 
GeV, and Gf = 1.16637 x 10"^ GeV-^, one gets 

I n.|2 

BR(K+ n+iyiUi) = 0.72 x 10"^^-^. (5.3) 

V 

us 

The function Di is proportional to the mass correction function Z)j[^], which is a function 
of mass ratios Xj = m^^/M^ and i/i = mfjM'y^. Neglecting the lepton mass in the value 
Dc for nic = 1.5 GeV is about 4 x 10^^. This value decreases with increasing lepton mass: 
for the r lepton with = 1.784 GeV, is reduced to 3.2 x 10~^. However, the value of 
Df increases almost linearly with rrit. In the range 80 GeV < rrit < 150 GeV, Dt is in the 
range 1.38 < Dt < 2.72. Taking a lower limit, rrit > 90 GeV, one gains a factor of at least 
400 in the ratio 

> 400. (5.4) 
With neglect of the QCD corrections, (|5.3|) becomes 

BR{K^ ^ 7r+z/,z7,) = 0.72 x lO-^lKdH - D,{x,) + ,^i^A(xt)P, (5.5) 

where the approximation V*gVcd ~ —V*gVud has been used. Neglecting lepton masses, one 
obtains the function D(x) as 



The weak contribution to the similar process Kl — ^ /i/i is given by 

BR{Kl ^ f^fiUak ^ 7.7 X 10-'>\Vus\-^\ReJ2V*VdC,{x,)\'', (5.7) 

i 

where Ci is defined by 

^/ \ 3 X 3x In X , , 

C{x) =x + + — -. 5.8 

4:1 — X 4(1 — xy 

As we have discussed in Chapter 4, the dominant contribution to the Kl — /i/i decay comes 
from the two-photon intermediate state. Taking the absorbtive part of the amplitude as 
calculated in ChPT, eq. ( [4.9| ), and assuming no interference between the real part of the 
same amplitude and the weak contribution as given in (|5.71 ), one obtains 



BR{Kl ixfi)disp ^{2±2)x 10-9. (5.9) 

The above assumption then yields the constraint 

\Re VtlVdCtixt)]^ < 1.7 X 10-3. (5.10) 

The constraint in ( p.lO| ) can be used keeping in mind that the possible interference in the 
dispersive part of the amplitude would weaken it. 
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5.3 QCD Corrections to the Kl fifi and tt^uu decays 

QCD corrections to the process tt+z/z/ have been discussed in a number of papersj 



13]' ^ controversy existed for some time, and has been resolved by Novikov et al. 

in [0. 

As far as the box diagram is concerned, there are three types of gluon contributions 
(diagrams (a) in Fig. 7): (i) gluons might be exchanged between d and s quarks, (ii) 
there is a quark-quark-W vertex correction, and (iii) there is a self-energy graph (mass 
renormalization) to an intermediate quark. 

It turns out that, working out QCD corrections in the Landau gauge, only the self-energy 
graph in Fig. 7 gives a logarithmic factor. 

The loop integral in the box diagram is of the form 



The last approximation is valid only in the LLA. 

An important point made by Novikov et al. is that the lower limit in ( p. Ill ) is not /i^ 



(arbitrary), but m^. For /i^ < m^, there are no logarithmic gluon corrections to the mass 
operator. This also solves the controversy with respect to previous calculations PU|- PT| . 

As is obvious from the graphs (a) in Fig. 7, the QCD corrections to the box graphs for 
Kl —>■ fifi and n^uu are the same. This is not the case for triangle graph corrections, 

to which we turn next. 

The triangle graph includes the T product of three hadronic currents. In contrast to the 
box graph, there exist now two independent distances. The calculation is simplified because 
there exists the Ward identity for the sdZ vertex. Therefore, one can reduce the calculation 
of the matrix elements of Kl — > fifi to the T product of two currents and the matrix element 
of the T product of two currents and the pseudoscalar density C75C. Finally, one gets that 
the bare result for Kl /U/i is multiplied by the Wilson coefficient 

+ 2^:"''" + nr"*"'), (5.12) 



cJ 



where ki = a<j(m^)/Q;s(M^). The expression (|5.12| ) is given for two families. Its gener- 



alization to three families is straightforward [p9[| . The following comments are in order: If 
as/47r ln(M^/m^) ~ 1, then the second term in ( 5.12 ) is of the same order as the neglected 



contributions in the first term. In the SVZ calculation [^] the second term has been kept for 
the following reasons: (i) its extrapolation to the free-quark limit is smooth, (ii) the leading 
term is not so large numerically since in the free-quark limit it vanishes, and (iii) this remains 
true for the three families as long as the t-quark is reasonably smaller than Mw- 

Generalization to three families is straightforward. 

The corresponding c-quark contribution is given by 

1 47r , 3 



""^ 2a,(m2)^ll 



-^-{kI'I^" - + ^.-^T" + 2Kr'^'' - 3«:r^^/^^). (5.13) 
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QCD Corrections for nit > 



The QCD corrections for the case rrit > My^/ are difficult to handle, since each of the 
sub leading mass terms gets differently 'renormalized'. However, the leading term for large 
mt is powerlike, i.e., the Wilson coefficients are simply 



C^M = C.p = ^, (5.14) 



the dominant contribution coming from the sdZ graph. This is , however, true for very large 
m^; in practice, we use in the range 80 — 200 GeV . Therefore, for large mt it is better to 
take the bare function Dt for C*p and similarly for C^^. 
The final expression then becomes 

BR{K+ ^ vr+z/.i.,) = 0.72 x lO-^lK^H - Ct, + ^^Cp. (5.15) 

The quantity |V^s^dP is a very important factor entering the calculations of rare decays. 
It appears as a 'weight' to the t-quark loop contributions. Using the measured K-M matrix 
elements, one infers from unitarity that 0.9985 < \Vtb\ < 0.9993 and 

\Vts\ < 0.054, \Vtd\ < 0.024, (5.16) 

which gives the constraint 

\Vtytd\ < 0.0013. (5.17) 

It is important to note that any violation of the relations (|5.16|) would be a signal of new 
physics beyond the standard model. The derived constraint in (|5.17| ) will be used to put an 
upper bound on BR{K^ tt+z/z/). 

In Fig. 8 we present the branching ratio for the decay tt+z/z/ as a function of m^. 

Since we have used the upper limit ( ^.17|) on K-M matrix elements, the result shown can be 



considered as an upper limit. One sees that varying mt in the range 90 GeV < mt < 200 GeV 
gives the branching ratio in the range 

BR{K+ TT+z/z/) = (1 - 4) X 10"^°. (5.18) 

The last Particle Data limit is[0 

BR{K+ Ti+uu) < 3.4 X 10"^ (5.19) 

The following comments are in order: 

(i) If one includes only the c-quark contribution (without the QCD corrections), one gets 
BR{K^ — > TT+z/z/) = 1.1 X 10~^^, for mc = 1.5 GeV. Adding the QCD corrections lowers the 
branching ratio to 3.1 x 10~^^. Clearly, both numbers are below the sensitivity of the BNL 
experiment E787 (a recent preliminary result is BR < 5 x 10"^, the expected sensitivity is 
2 X 10"^° ). 
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(ii) For larger values of rrit, the t-quark contribution dominates and the results presented 
in Fig. 8 can be used in different scenarios. If the experiment measures considerably larger 
values than predicted, it is unhkely that it could be explained by pushing mj too high, since 
that would violate the Kl fift bound ( [4.3|) . However, it would be possible to violate the 
unitarity constraint in (|5.171 ) because, for example, of the existence of the fourth family. 
Therefore, any clear result beyond the predicted branching ratio should be considered as 
a serious signal of new physics. Clearly, further improvements of the upper limits on rrit 



and/or |Vj*Vtd| would further constrain the range (|5.18 



(iii) In the case that rrit is determined from experiment, the eq. (|5.15| ) may be used to 
test the standard-model prediction for the branching ratio. On the other hand, the measured 
branching ratio ( [4.3|) may be used to derive an upper limit on rrit. Both possibilities are, of 
course, valid provided no significant impact of 'new physics' is present. 

(iv) We have kept QCD corrections for the c-quark contribution only. Since the latter is 
not dominant, the errors entering through the LLA, the charmed quark mass, etc., are not 
significant for the branching ratio. For example, the threshold effects of quarks have recently 
been calculated for C^p and they slightly modify the resulting branching ratio. 

(v) Recently, a penguin-box expansion has been performed systematically [p^ ] ; cf. also 
[p5| . It enables one to perform a very detailed analysis of K ^ fifl and tt+z/z/ with 
constraints coming from the e, e' analysis included, yielding BR{K^ tc^uu to be around 
1 X 10"^°. 



5.4 TT^uiy decay in the minimal supergravity model 

We first supersymmetrize the box diagram, i.e., the internal quarks and W-bosons become 
squarks and gauginos (Fig. 9a). This has to be compared with the box diagram in the 
standard model which is given as 

A's^iK-^ - n+i^u) = -^^V*VtM^uyi)-9ix^,yi)]i-SL7Mii^ri^L), (5.20) 



where g{x,y) is given as 



gix,y) = J—{-(^)Hnx + i^f\ny+^ - -^}. (5.21) 

y — X X — I X — I y — i- 



The supersymmetrization leads to the following expression: 

AZy{K^ ^ ^^^^) ^ ^^^tsVtdWum) - giiuMisLlMiv^^VL). (5.22) 

w 

Numerically, the ratio 

'^=\Ausy/AlA (5.23) 

is of 0{1). We have used m| — m\ ~ and m^r = 40 GeV , mi = 40 GeV. If one fixes 
rriu = 70 GeV, then 3? grows with rrit; for rrit = 90 GeV, 3? = 0.4, and for rrit = 160 GeV, 3? 
reaches 1. For very large values of rut, e.g., rut = 200 GeV, ^ = 1.4. With increasing ruu, 3? 
becomes somewhat smaller: for rUu = 100 GeV, rut = 60 GeV, 3? = 0.2, and reaches 1 for 
Tilt = 200 GeV. 
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Although the SUSY contribution may, in principle, be a significant correction for large 



values of nit, it is unlikely, as noticed by Bertolini and Masiero|^^, that this really would 
happen. In fact, the large 3? would mean the large SUSY contribution to the /i/i box. 

This is even more pronounced since the Dirac algebra suppression in the fifi box by 

factor of 4 does not appear in the SUSY version and, morever, sneutrinos that enter the 
Kl —>■ box further increase the amplitude. 

The promising contributions appear to be spenguin graphs (b) and (c) of Fig. 9. Graphs 
(b) with the external photon contribute, for example, to the 7r+e^e~ decay. With the 

external Z-boson, graphs (b) and (c) contribute to the Kl /x/i and tx^vv decays. 

For set (b) it turns out that the self-energy graphs exactly cancel with the vertex correction 



graphs. Therefore, one has to include two L-R insertions in the squark propagator ||47||. The 
resulting mixing (11 — da is given by Am3/2Md (cf. ( ^.14D ): 



m 



'3/2 



The standard-model value for /i/i is given by 

G p 



(5.24) 



AsuiJ^L /i/i) 



a 



y/2 TT sin^ 6 



(5.25) 



with C{xt) given by (|5.8| ). The SUSY contributes as||47 



'72^ 



V*Vu{Ams/2m,/mlY[C{x,) - (^)'C(xrf)](si7t^L)(/iL7>L) 



rrib 



with C{x) defined as 
C{x) 

This leads to 



x-iy'2 x-1 



2«., 



(x-1] 



In X — 



1 



^{Kl /i/i) < 77— sin 6^{Ams/: 
6 a 



x{x — 1 

mb.2, C{xb) , 
'■ml' ^CixtV 



-}• 



(5.26) 
(5.27) 

(5.28) 



The corresponding expressions for A{K~^ vr+z/z/) can be obtained by changing C{xt) to 
D{xt), i.e. 

sR(i^+ ^ ^+^^) = \^^mKL ^ /i/i). (5.29) 

Numerically, (|5.29|) presents a sizeable correction only for small values of the supergravity 
parameters. For 7/13/2 = 100 GeV , nig = 70 GeV, mi = 40 GeV, and choosing mj = 90 GeV, 



0.3. 3? decreases with increasing m^; for mt = 150 GeV, 9ft = 0.14. On 



the other hand, larger values of suppress 3ft significantly: for m^ = 60 GeV, mt = GOGeV, 
3ft is lowered to 0.1. 

Comparison with the supersymmetrized box shows that for mt = 90 GeV, both super- 
symmetric contributions are of the same order, 30 — 40%. With increasing mt, ^^p<^^9uin 
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diminishes and 3?''°^ increases. Although a precise prediction is difficult because of the arbi- 
trariness of SUSY parameters, it seems that in the mass range 90 GeV < mt < 150 GeV the 
SUSY effects^ may reach 50 — 100%, leading to the enhancement in the branching ratio: 

BR{SM + SUSY) ~ (2 - A)BR{SM). (5.30) 

Prospects to detect the supersymmetry in this rare decay are presently not very bright. 
However, once and/or K-M matrix elements are determined, the prediction ( |5.15| ) for 



the branching ratio in the SM is rather reliable. The precise measurements may be able to 
disentangle the possible departures from the standard model. For example, the branching 
ratio enhanced by a factor of 4 may be due to supersymmetry, but also to the fourth gener- 
ation; the precise determination, however, of the K-M matrix elements plus unitarity could 
eliminate the latter possibility. 

5.5 7r+77 decay 

The i^'"'"-decay into a pion and a pair of the lightest superparticles is a very interesting 
decay mode since it belongs to the 'direct'-decay type: the superparticles are produced as 
real particles. However, in the minimal = 1 supergravity model we are discussing, it 
is not likely that this mode will appear for the simple reason: this decay is not allowed 
kinematically. The underlying assumptions are: that the 7 is the lightest SUSY particle 



and i?-parity is conserved. The existing limits are the following ||19||: m^y > 15 GeV for 
mj = 100 GeV, > 5 GeV for mg = 55 GeV, and there is no lower limit for mg > 58 GeV. 
Although for the heavy selectron there is no lower limit on m^, there is a cosmo logical 



limit []T^ that requires the photino to be heavier than 300 — 500 MeV or very light, m^y < 



100 eV. Actually, there is a new allowed region m^y = 4 — 20 MeV provided the bound 
"^3/2 < 40 TeV is satisfied. 

The second underlying assumption (R-parity) may also be relaxed. The SUSY models 
with i?-breaking (which allows the lightest superpartner to decay) have been considered for 



a long time and recently discussed again Therefore, it looks meaningful to discuss the 
possibility of the 7r"'"77 decay both from the experimental and theoretical point of 

view. 

The largest contribution to the decay comes from the presence of large effects of flavor- 
changing quark-squark-photino couplings, which simply leads to a tree-graph s^di —>■ 77 via 
a c?^ exchange. The result is very simple P7[] 



A{K+ ^ 7r+77) = -^[cm?/(m^/2 + cmX/2)]V^i:V^td(^L7A)(77"757), (5-31) 
and has to be compared with the dominant part of the n^uu decay: 

A{K+ ^ 7r+z/z7) = £—^V*VuDt{xt){sL'yMi^LY''L)- (5.32) 

V2 vr sm Oy, 



^^There is a contribution of the graph where a nutral higgsino plays the role of giuino; this contribution 
is, however, suppressed by small Yukawa couplings. This is not the case for charged higgsinos: although 
enhanced, this contribution is still smaller than the contribution of the standard model[[t7|. 
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It is easily seen, that the former decay may well dominate over the lattei|3 

Therefore, the experimental upper limit on tt^+ ^nothing^ may be converted to 

the lower limit on the squark mass which appears as m^^ in the expresion ( |5.31D . Using 
the present experimental bound in (|5.19|) one derives the following constraint on the squark 
mass: 

rrig > 50 GeV. (5.33) 

As we have said, this limit is valid provided, that the decay — > n'^'y'y is allowed kine- 
matically^. 



13 



In order to compare both amplitudes, ( 5.32| ) has to be muhiphed by a/3 for 3 neutrino species (the 
Majorana nature of the photino is already taken into account in (5.31) ). 

The tree-graph decay into neutral higgsinos is suppressed for the same reason as in the spenguin 
contribution involving higgsinos. 
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6 Constraints on Supergravity Parameters from CP 
Violation and K ^ tttt Decays 

6.1 - Mixing 

In this section we discuss how the Kl — Ks mass difference and the e parameter in the 
— system constrain a general version of the minimal = 1 supergravity model 142 



The supersymmetric contributions to Aitlk and e are box graphs depicted in Figs. 10 and 
11. The graphs in Fig. 10 represent the mere supersymmetrization of the 'standard box': the 
internal quarks are replaced by their scalar partners, the squarks Wj, whereas the W bosons 
and the Higgs scalars are substituted by their fermionic partners, W and if, respectively. 
The gluon exchange diagrams in Fig. 11, proportional to the strong coupling, are expected to 



give the dominant contribution [g^, 0. One usually assumes no great cancellation between 
different contributions to Am^^ and e, and requires each contribution to be less than the 
experimental values[^ 

Diagrams (a) in Fig. 10 yield a constraint [ 42 

(V^^V),, < (6.1) 

m3/2 100 Mw 

which implies near degeneracy amongst the squarks. This is consistent with the usual as- 
sumption about masses of scalar fields of chiral superfields 

M!, = m\i^b,, (6.2) 

at a renormalization point of the Planck mass. 

The contributions of the diagrams (b) in Fig. 10 place a constraint on the flavor matrices 
iu, and io 

iu,D = m3/2AXu,D + iu,D (6.3) 

appearing in the trilinear scalar couplings (|2.5| ). The ^u^d is an arbitrary small matrix. 
However, more stringent constraints on ^ujj are set by diagrams (a) in Fig. 11, leading 
to effective operators with mixed L-R helicities, i.e., with the structure typical of penguin 
graphs in the SM. The local AS" = 2 hamiltonian generated by this graph is given by 

neff = -^a^s^^H^fi^m - 3^2), (6.4) 
120 m°/2 "^3/2 



where 



and 




sll.di -s'kYd'n I I (6.5) 



-3 



/(x) = 20 J^C C'(l-C)[C + a:(l-C)]" 

/(I) = 1. (6.6) 



Franco and Mangano|Q gave the general constraint assuming a coherent contribution of diagrams (a) 
and (b) in Fig. 11. 
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The notation for the mass parameters in ( |6.4| ) is explained in . 



In order to estimate the supergravity contribution in Aitlk and e, we need an estimate of 
the matrix elements of the operators Oi and 6*2- Using the vacuum saturation as a starting 



point, we define [46 



< Ky,\K' >= \fWA^ + ^)Be.. (6.7) 



" c 



With this definition, the vacuum saturation values are Bq^ = —9.6 and Bg^ = —4. As usual, 
the L-R helicity structure of our operators disqualifies the reliability of the vacuum saturation 
estimate. A better estimate is obtained by making use of the QCD duality approach |56[ . 
We have to study the behavior of the two-point function 

^g^{q^) = I J d^xe'"-'' < O\T{02{x)el{O))\O > . (6.8) 

The first step is to calculate ^Im '4'92{t) both in QCD using the representation ( |6.5| ) of 6*2 
and in ChPT using the chiral representation 

QcMrai ^ Be,^{&) : {^S.U^,U%s {^^-fU^^^UU ■ ■ (6.9) 

Next, we establish the duality region in terms of FESR's. One needs two sum rules to 
fix So, which is the onset of the asymptotic QCD continuum. These are 



rso 1 

3^0 = / dt—Imipg^it), 



'■K 

= / dtt-Im^e^{t). (6.10) 

Once So is fixed, any of the above sum rules leads to the value of Bg,^ . The ratio of the sum 
rules 

^ = ^^ 6.11 

5so Mq 

is a very sensitive test of duality. It does not depend on Bg^ and may be used to fix so. 

Fig. 12 shows the results plotted in both in ChPT and in QCD. The departure from 
the asymptotic QCD prediction (the dashed line) is due to mass and condensate corrections. 
The dots are the values of r in ChPT, and the solid line is the full QCD result. The duality 
region is clearly established in the range 8 < so < 11.5 GeV^. Any of the sum rules in (|6.10| ) 
leads to the value of Bg^ once sq has been determined. 

Fig. 13 shows the results obtained from the sum rule 3?o- Bg^ plotted as a function of Sq 
shows a plateau behavior in the duality region. The results obtained from the sum rule 3?i 
are similar. The values for Bg^ agree within 2% in both cases. Taking into account various 
theoretical uncertainties [HBIl, we obtain the final estimate 



1 1 = 0.25 ±0.15. (6.12) 

There are some principal difficulties in determining the matrix element of the operator 
= 6i — 362- The QCD duality approach determines only the absolute value of the B 
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parameter and a separate determination of Bq-^ would not be very useful. However, the 
operator 6i is of higher order in ChPT with respect to 62 and is expected to be significantly 
smaller. In addition, the QCD spectral functions of the operators 6 and 302 show very similar 
t-dependence and are of about the same magnitude ||46||. Therefore, we may neglect the 61 
contribution to a very good approximation. 

Using our estimate (|6.12| ) of the matrix element of the 62 operator and the experimental 



values for Attik and e, we obtain the following constraints on and M'^, 



DS 



M^r^M*f^^, .ml/2 
Re{ <4.5x IQ-^-ff,, , ^ 



A/r2 i^(„„2 )] ' 
2 _2 



Ml^M*f^^, 7^^3/2 r./ 



m 



3/2 



m 



(6.13) 
(6.14) 



3/2 



These constraints are rather sensitive to the variation of the ratio Tng/m^i2- If one takes 



rrig = m3/2, one finds that ( 6.13 ) and ( 6.14 ) are weaker by a factor of 5 than the corresponding 



constraints of Dugan et al. 



The superbox diagrams (b) in Fig. 11 yields the effective hamiltonian ||35|| , - pO 



(6.15) 



where the sum runs over indices u, c, and t, and 



m 



Xi 



f, X2 



m 



2 ' 



3^3 



b_ 

ml 



(6.16) 



2 

mi 



m- + mj. The function S{x,y) is 



The down squark masses are related to each other through the equation ( p.l7|) . Thus, one 
may assume that in a good approximation m| = 
given by 

^^ Kix)-Kiy) 



S{x,y) 



X 



y 



X 



(6.17) 



K{x) 
Kx) 



x^lna; 1 
+ 



(l-x)2 
X Inx 



1 — a; 
1 



(1 — xy 1 — X 

Using the unitarity of K-M matrix one finds the super-box contribution to Amx 



(6.18) 



(Am 



K ) sbox 



54m? 



BfKmK\V*Vtd\'AS, 



where 



AS* = S{X3,X3) + S{Xi,Xi) - 2S{Xi,X3), 



(6.19) 



(6.20) 



and B = 1 for the vacuum saturation estimate of the LL operator matrix element between 
the and K'^ states. For m^ <^ mg, AS can be replaced by the second derivative of S{x, y) 
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and ( p.lTp simplifies, yielding an mf dependence of Aitlk and e We shall, however 



use the exact form (|6.19|) . 

Eq. ( |6.19| ) can be used to place lower bounds on the gluino mass and down-squark masses, 
provided one has sufficient information about the quark masses and K-M parameters. Un- 
fortunately, the top-quark mass is still not known. The best one can do is to use the 



present experimental lower bound rrit > 91 GeV|T^ and the upper limit rrit < 180 GeV 



from the standard model constraints |^6[. As regards the K-M parameters, a rigorous upper 
bound given by ( |5.17] ) can be used. The lower bounds on rrig and for two values of mt 
are as follows 



mt 100 150 

rrig > 40 50 

> 105 160 

The imaginary part of the strong superbox diagram (b) in Fig. 11 yields a contribution 
to the e parameter. Whereas for typical values nig ~ ~ "^3/2 = 100 GeV {AmK)sbox is 
rather small (of the order of 0.1(Amx)exp ), the contribution to e can be quite large p5|. By 
making use of ( |6.15| ), one finds that 

I I 2BfK^K 2 2/ NAC (an^\ 

\eshox\ = TT^Oi^^^ S^S2[S2SzSs)/\b, (6.21) 

iUo TJlglATrij^ 

which compared with eexp = 2.27 x 10~'^, and assuming {AmK)sbox/{A.mK)exp ~ 0.1, leads 
to a stringent constraint on S2S3S5: 

S2S3S5 < 2 X 10-^^. (6.22) 

•^2 

This bound has important consequences for the contribution of penguins and superpenguins 
to e', which we discuss in the next section. 

6.2 CP Violation in the = 1 transition 

Supersymmetric contributions to e' are represented by the superpenguin operators depicted 
in Fig. 14. Again, the corresponding diagrams involving winos are neglected because they 
are proportional to rather than to g"^. The superpenguin (a) in Fig. 14 leads to a phase 
in the K]^ — > tttt amplitude, ^spen compared with ^pen from the standard model equals P2| 



^spen 1/5^5x2 {T^t/m^ 



\2 



dJ 



ipen 5 g ln(mt/mc)^ 



(6.23) 



This could be larger than 1 if nit is very large. However, in this case, the superbox diagram 
dominates e and yields a strong bound on 5253^5, ( |6.22| ), which in its turn makes the penguin 
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contribution very small. Thus, the conclusion is that the superpenguin (a) in Fig. 12 may 
give a significant contribution to e'/e only in a very small region of parameter space|4^ . 



The scalar mass insertions and which have been discussed in the preceding 

section, also appear in the AS" = 1 transition dipole moment operators]^, represented by 



the superpenguin diagram (b) in Fig. 14. The corresponding effective hamiltonian is given 
by 

'Heff = i^^J'{x)[^sna,,F^''dL + ^sa^^F^'^dni (6.24) 
327rm3/2 "^3/2 "^3/2 



with 



"^3/2 

T{x) = ^[l2^{x) + \-^%-)l ^(1) = 1 + 1, 

%x) = f dai - C?\C + x{l - C)]-^ (6.25) 

This effective hamiltonian could give a significant contribution to e' and to the A/ = 1/2 
amplitude in the K ^ -k-k decay. Using, for example, the bag model estimate[^ of the above 
transition moment operators one finds the phase and the magnitude of the A J = 1/2 
amplitude 

itransmom = 250 ):F{x)Im + -^), 6.26 

m3/2 "^3/2 

ailns^om = 1-2 X 10-5 )J^{x)Re + -^). 6.27 

m3/2 "^3/2 



The constraints ( 6.26 ) and ( 6.27 ) allow us to assume that the imaginary parts of Mlj^ and 



Mq^ are much smaller than the real parts. Assuming further that Re M^d ^ Re Mf)^, one 
obtains 

\e/e\trans mom < 2 X 10-'<l>(a;), (6.28) 
allLmom < 20 X 10-^^x), (6.29) 

with 

$(x)=^(x)[/(x)]-i/2 (6.30) 



In view of the recently measured values i?e(e'/e) = (2.2 ± 1.2) x 10 ^ [|1^], the bound (|6.28| ) 



is not satisfactory. On the other hand, the bound (|6.29|) , compared with the experimental 
value a^/^ = 27 x 10^^ GeV, opens an interesting, although perhaps unreahstic, speculation 
that a great deal of yet unexplained A J = 1/2 enhancement could be attributed to the 
transition dipole effective operators induced by the extended supergravity model. 

The bag model estimate may be rather crude for the operators containing gluon fields. The QCD duality 
approach is, however, very difficult to apply here because of the huge corrections, typical of AI = 1/2 
processes] 62 . 
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The mass ratio mg/m^i2 is yet unknown. Assuming = 7713/2 and requiring 

1^ / ^\trans mom ^^^trans mom ^ 3 X 10 . (6.31) 

one finds the following constraint 



,„(iHk + £|£)<8xlO-'(-^). (6.32) 
This constraint is more stringent than the one obtained from ( |6.26|) 
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7 Conclusions 



The processes discussed in this review have been selected in such a way that they are inter- 
esting from both the theoretical and experimental point of view. The latter means that rare 
processes we are discussing are expected to be detected in the near future. Of particular 
theoretical interest are processes which either have significant corrections (here 'significant' 
means corrections of the order C(l) ) owing to the supersymmetry, or the standard model 
predictions are clear and unambiguous. In the latter case, the departure may indicate physics 
beyond the standard model, and in the former case precise measurements again check the 
possible new physics. 

The process — > ti^vii receives supersymmetric corrections which are of the order 1. 
The resulting branching ratio is larger by a factor of 2 — 4 with respect to the standard model 
prediction. The latter is a function of the K-M parameters V^gVtd (only the upper hmit is 
known), and of the top-quark mass, for which we know the lower limit. Partial knowledge of 
the above parameters leads to the prediction BR{K^ tt^vv) = (1 — 4) x 10~^°. The range 
given will become more narrow if the lower limit on mt increases. This decay is dominated by 
short-distance dynamics, and long-distance effects are expected to be small, and its standard 
model prediction is relatively free of usual uncertainties. Therefore, if the decay is measured, 
it can also be used as a way to determine mt- 

Direct production of superparticles (photinos) is possible in the decay 71+77, 
provided it would not appear to be kinematically forbidden. The contribution is large enough, 
so that the present upper bound on K'^ — > 7T+ 'nothing^ leads to the constraint mj > 
50 GeV. 

Supersymmetric contributions to ArriK and e may be significant for large values of mt- 
This yields constraints on gluino and squark masses as well as the L-R helicity mixing 
mass parameters M^g and M§d. Further constraints on these parameters follow from the 
requirement that the contribution to e' of the transition dipole moment operators in the 
general version of the supergravity model must not be too large. In the constrained version 
of the model the superpenguin contribution to e'/e appears to be rather small. 

Given the fact that supersymmetric effects are rather tiny, it is important to reduce 
uncertainties in the calculation, especially the ones coming from the nonperturbative (con- 
finement) effects in QCD. Some recent approaches, such as chiral perturbation theory, the 
large-A^c expansion, QCD hadronic duality sum rules, lattice QCD, are discussed. Partic- 
ularly interesting appear predictions of ChPT in rare decays. Especially, the reachness of 
radiative K decays clearly shows the predictive power of ChPT. Some of these predictions 
are to be tested in the forthcoming experiments in the near future. 
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Figure Captions 



Figure 1. Typical graphs in ChPT. (I) Tree-level decay, (II) pion decay constant and wave 

function rcnormalization, (III) loop corrections to the — mixing, (IV) loop corrections to the 
decay. The square is an insertion of the weak lagrangian and the circle is a strong interaction 
vertex. 

Figure 2. The 'eight '(a) and the 'eye'(b) diagram contribution to the correlation between two 
mesons and a four quark operator. 

Figure 3. Rare-decays graphs at the quark level, (a) box graph for the Kl /U/v, decay, (b) 
box graph for the Kl — > 77 decay, (c) sdZ-vertex graph for the Kl fijl decay, (d) two-photon 
contribution to the Kl — fJLji decay, (e) box graph for the K^ — ir'^uu decay. 

Figure 4. Rare decays in ChPT. (a) two-photon contribution to the Kl — ^ jJ-Jl decay, (b) loop 
graphs for the Kg — > 77 decay, (c) pole contribution to the Kl 77 decay. 

Figure 5. Supersymmetric contributions to the /x — *■ 67 and g — 2. (a) /x — *■ 67 decay in SUSY, 
(a), (b), and (c): contribution to g — 2. 

Figure 6. K~^ n'^uv decay, (a) box graphs, (b) sdZ-vertex contributions. 

Figure 7. QCD corrections to the K^ — > tt^uU decay, (a) corrections to box graphs, (b) 
corrections to the sdZ-vertex graph. 

Figure 8. Upper limit on BR{K^ tt'^uU) in the standard model as a function of top-quark 
mass. 

Figure 9. Supersymmetric contributions to the K~^ tt^vu decay, (a) supersymmctrizcd box 
graph, (b) supersymmetrized sd'y and sdZ vertices, (c) super symmetrized sdZ vertex with mass 
insertions. 

Figure 10. Supersymmetric box graphs proportional to weak couplings. 
Figure 11. Supersymmetric box graphs proportional to strong couplings. 
Figure 12. The ratio of the sum rules r plotted versus sq in ChPT (dots) and in full QCD 
(solid line). 

Figure 13. Results for Bq^ as a function of sq. 
Figure 14. Supersymmetric penguin graphs. 
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